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We present a search for exclusive Z boson production in proton-antiproton collisions at ^/s = 
1.96 TeV, using the CDF II detector at Fermilab. We observe no exclusive Z £'^£~ candidates 
and place the first upper limit on the exclusive Z cross section in hadron collisions, crexcl{Z) < 
0.96 pb at 95% confidence level. In addition, we observe eight candidate exclusive dilepton events 
from the quantum electrodynamic process pp p'yyp pt^l~p, and measure the cross section for 
Mu > 40 GeV/c^ and |r/^| < 4 to be cr = 0.24to io pb, which is the first measurement for this mass 
range and is consistent with the standard model prediction. 
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PACS numbers: 



At the Tevatron pp collider it is possible to produce Z 
bosons exclusively, in association with no other particles 
except the p and p: pp pZp. The colliding hadrons 
emerge intact with small transverse momenta, pt [H- 
The process is predicted by the standard model (SM) to 
proceed via photoproduction. A radiated virtual photon 
fluctuates to a qq loop which scatters elastically by two- 
gluon exchange on the (anti)proton and materializes as 
a Z^ as shown in Figure |l(a)[ The same mechanism gives 




P P 



P P 




(a) 



FIG. 1; (a) Exclusive photoproduction of a, Z boson and (b) 
exclusive dilepton production via two-photon exchange at the 
Tevatron. 

photoproduction of the vector mesons J/tp,ip{2S) and 
T, which have been studied in ep collisions at HERA ^] 
and recently observed for the first time in pp collisions 
by CDF Q . The SM cross section for exclusive Z pro- 
duction is predicted to be aexci{Z) = 0.3 fb ^Sj, and is 
thus below the threshold for detection. An observation 
at the Tevatron would therefore be evidence for beyond 
SM (BSM) physics. A BSM theory of the pomeron P [g] 
predicts a much larger cross section, possibly orders of 
magnitude larger, but without a quantitative estimate. 
In this theory the pomeron couples strongly to the elec- 
troweak sector through a pair of color sextet quarks which 



state particles in exclusive dilepton events are identical 
to those in exclusive Z production with leptonic decay, 
the only difference in the signature being the Mu distri- 
bution and other kinematics. Agreement with the precise 
theoretical prediction therefore gives us confidence in our 
sensitivity to selecting exclusive Z bosons. 

CDF II is a generalpurpose detector which is described 
in detail elsewhere [8|; here we give a brief summary. 
Surrounding the collision region is a tracking system 
consisting of silicon microstrips and a cylindrical drift 
chamber, the central outer tracker (COT), in a 1.4 Tesla 
solenoid. The tracking system tracks particles with 
Pt ^ 0.3 GeV/c and pseudorapidity |?7| < 2 Central 
and end-plug calorimeters cover the range \r]\ < 1.3 and 
1.3 < |?/| < 3.6 respectively, with separate electromag- 
netic (EM) and hadronic (HAD) compartments. Outside 
the calorimeters, drift chambers measure muons in the re- 
gion < 1.0. The regions 3.6 < |?7| < 5.2 are covered 
by lead-liquid scintillator calorimeters called the mini- 
plugs Q. At higher pseudorapidities, 5.4 < \r]\ < 7.4, 
scintillation counters called beam shower counters (BSC) 
are located along the beam pipe. Gas Cerenkov detec- 
tors covering 3.7 < \ri\ < 4.7 measure the luminosity with 
a 6% uncertainty by counting inelastic interactions [lo| . 
Tracking detectors in a Roman pot spectrometer 11| , lo- 
cated 57 meters from the interaction point, can detect 
antiprotons with small and 0.03 <, ^{p) <, 0.08, where 
^(p) is the fractional momentum loss of the antiproton [3]. 
These detectors were operational for approximately 30% 
of the data set used in this analysis. 



For the 77 — > event selection we select a sample 

of £'^£~ pairs in a kinematic region where this process has 
not previously been observed, with Mu > 40 GeV/c^ and 
lepton transverse momenta pfp > 20 GeV/c. For the ex- 
clusive Z search we select a subsample with an invariant 
This Letter presents a search for exclusive ' Z produc- ^^^^ ^^'^^'^ *o ^ ^'^ss, 82 < M« < 98 GeV/c^, and 



contribute to the quark loop shown in Figure 1(a) 



tion with the Z decaying to a /i /i~ or e"'"e~ pair, 
measurement of the cross section for exclusive /i"*"^' 
e'^e~ production with dilepton invariant mass Mu 



and a > 25 GcV/c. The /i^/i events are collected with a 



and 
> 40 

GeV/c^ and \Tji\ < 4. We use the CDF II detector at 
the Tevatron with pp collisions at a center of mass en- 
ergy y/s = 1.96 TeV. The exclusive dilepton process is 
expected in quantum electrodynamics (QED) through 

For 



P77P — > p£^£ p, as shown in Figure 1(b) 



PP 

the remainder of this Letter we will refer to this process 
as 77 £^£^ for convenience. We have previously ob- 
served 77 e+e~ with 10 < Mee < 40 GeV/c^ Q and 
77 fJ-^fJ-'' with 3 < M^^ < 4 GeV/c^ [1] and measured 
cross sections in agreement with expectations. The final 



'Deceased 



trigger requiring one muon with px > 18 GeV/c. Offline 
we require two candidate muons. One muon must be de- 
tected in the COT, the central calorimeter, and the muon 
chambers, and therefore has jry^j < 1.0. To increase ac- 
ceptance the second muon is only required to be detected 
in the COT and therefore has \r]fj,\ < 1.5. Events consis- 
tent with cosmic rays are eliminated with an identifica- 
tion algorithm J^] that uses the timing of the COT drift 
chamber hits. The muon kinematics are found from the 
COT track momentum measurement. The e~^e~ events 
are collected with a trigger requiring one central electron 
with Pt > 18 GeV/c. Offline we require one candidate 
electron reconstructed in the central EM calorimeter and 
matched to a COT track, and a second electron recon- 
structed either in the same way or in the end-plug EM 
calorimeter where no matching COT track is required. 
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since the tracking efficiency is lower in this region. The 
central electrons have |?7e| < 1.3 and the end-plug elec- 
trons have 1.3 < |77e| < 3.6. The electron kinematics are 
found from the calorimeter energy measurement, but if a 
track is matched to the calorimeter cluster it is used to 
determine the electron direction. If no track is matched 
the z position of the interaction is measured from the 
other electron track, and is used to determine the kine- 
matics. Events with two central electrons are denoted 
CC e+e~ events, and those with one central and one end- 
plug electron are denoted CP e^e~ events. The \jr , 
CC e+e~, and CP e+e" events are each treated as sepa- 
rate final states which are ultimately combined together. 
With an integrated luminosity of 2.20(2.03) fb^^ in the 
electron(muon) channels we find a total of 317,712 candi- 
date dileptons with M« > 40 GeV/c^, of which 183,332 
are in the Z region 82 < M« < 98 GeV/c^. 

Starting with the dilepton samples we select events 
that are consistent with arising from exclusive produc- 
tion, by requiring that no other particles are produced in 
the collision. We veto events where any additional tracks 
are reconstructed in the COT or the silicon tracker, or 
where any of the calorimeters have a total energy depo- 
sition above that expected from noise. For this purpose 
the calorimeters are divided into five sub-detectors and 
the energy of all towers is summed, excluding those tra- 
versed by and surrounding the charged leptons, to give 
five TjE values. Each YjE is required to be less than 
a threshold, which is determined by studying two con- 
trol samples: (1) events selected with a random bunch- 
crossing (zero bias) trigger with no tracks in the event, 
which should give distributions dominated by noise and 
(2) W ^ iiy events with no detected tracks other than 
one from the charged lepton, which should give the dis- 
tributions expected for non-exclusive Z £~^£~ events 
with no additional tracks. The production mechanism 
for non-exclusive W bosons is very similar to that for Z 
bosons and the cross section for exclusive W production 
{pp — > n Wp) is negligible, making them an excellent con- 
trol sample. The chosen energy cuts are Y,E < 3 GeV in 
each of the East and West miniplugs, < 5 GeV in both 
EM plugs, < 7 GeV in both HAD plugs, and < 0.35 GeV 
in the central EM calorimeter. 

These exclusivity cuts reject exclusive events that are 
in coincidence with additional inelastic pp collisions. It is 
therefore necessary to define an effective integrated lumi- 
nosity J Ceff- The fraction of bunch crossings, selected 
from the zero bias trigger, that pass the exclusivity cuts 
is used to establish that / Ceff is 20.6% of the total in- 
tegrated luminosity. The efficiency is found from distri- 
butions reweighted to account for the difference in the 
instantaneous luminosity profiles between the zero bias 
events and the Z events. As a cross check we also estimate 
/ Cef f using Poisson statistics and the mean number of 
expected interactions per bunch crossing as a function of 
instantaneous luminosity to be 18.7% of the total inte- 



grated luminosity. Since the method using the zero bias 
data properly takes into account events with no interac- 
tions that fail the cuts due to noise in the calorimeters 
and fake reconstructed tracks, and events with a very soft 
interaction that pass the exclusivity cuts, we use 20.6% 
to determine / Cef f , and take the 9% difference between 
the two methods as a systematic uncertainty. We find 
J Ceff = (403 ± 45) pb-i and (467 ± 50) pb"! for the 
fj-^ fj.~ and e+e~ samples, respectively, where the uncer- 
tainty includes contributions from the difference between 
the two methods and the 6% uncertainty on the CDF lu- 
minosity measurement. 

In order to reduce the background from 77 — > £~^£~ 
events where the proton dissociates into forward-going 
hadrons, we also make cuts on hits in the BSC detectors. 
An event is vetoed if any photomultiplier has hits above 
threshold, which occurs in 76% of zero bias events that 
pass all the other exclusivity cuts. This inefficiency is 
included in the acceptance. 

A total of eight events pass the 77 t^£^ selec- 
tion criteria and no events pass the tighter exclusive 
Z — > criteria. We use these events to measure the 

cross section for the 77 £~^£~ process and we set an 
upper limit on the cross section for exclusive Z produc- 
tion. To do this we need to determine the acceptance 
for reconstructing the events and the expected number 
of background events. 

We calculate the acceptance for reconstructing 77 
e+£~ events using the LPAIR [i3| Monte Carlo (MC) 
event generator together with a GEANT IJ] simulation of 
the CDF detector. We apply corrections to account for 
changes in the acceptance due to internal Bremsstrahlung 
from the leptons, using the PHOTOS [l^l MC event gener- 
ator. The acceptance for the exclusive Z search is found 
from the pythia [3| MC event generator, which simu- 
lates non-exclusive Z/j* £^£^ events. Corrections are 
applied to account for the difference in kinematics be- 
tween non-exclusive and exclusive production. We con- 
sider the Z pt distribution, which is assumed to be be- 
tween and 2 GeV/c for exclusive Z production, the Z 
rapidity yz distribution, obtained from Ref. and the 
angular distribution of the leptons. The latter is assumed 
to be that for the decay of a spin-1 boson into two spin- 
1/2 fermions, which has the form (1 -I- cos'^ 6**), where 9* 
is the polar angle between the outgoing lepton and the 
proton direction in the boson rest frame. 

The backgrounds to 77 — > £~^£~ events are non- 
exclusive Z/7* £'^£^ events that pass the exclusivity 
cuts, and 77 — > £~^£~ events where the proton or an- 
tiproton dissociates and the products are not detected 
in the forward detectors. The former is found to be 
0.28 ± 0.19 events by assuming the fraction of non- 
exclusive Z/j* — !■ £~^£~ events passing the exclusivity 
cuts to be the same as that for non-exclusive W ^ £v 
events. This fraction is found from W ^ £v data sam- 
ples, selected by requiring a high pT lepton and large 
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missing transverse energy, to be (9 ± 6) x 10~^, where 
the uncertainty comes from the statistics of the sam- 
ples. The latter is found from the lpair event generator, 
which also simulates 77 — > events where one or both 
(anti)proton dissociate. We use the minimum bias Rock- 
efeller MC [l3|, which fragments the excited (anti)proton 
into a nucleon and pions, to predict the fraction of dis- 
sociation events that fail our exclusivity cuts due to par- 
ticles in the region \ri\ < 7.4, which is the edge of the 
BSC acceptance. We predict a total background of 1.45 
± 0.61 events, where the uncertainty comes from vary- 
ing the exclusivity cuts and observing how the number 
of events changes. 

The backgrounds to exclusive Z events are non- 
exclusive Z/j* — > £^£^ events that pass the exclusivity 
cuts and exclusive 77 — > £'^i~ events that have Mu in 
the Z mass window. The former is found to be 0.163 ± 
0.099 events using the method described above, and the 
latter is found from the lpair MC samples to be 0.492 
± 0.061 events. We do not include a dissociation back- 
ground for the exclusive Z search; instead we quote an 
upper limit on the cross section for a Z produced with no 
other particles with jryj < 7.4. 

There is no evidence of cosmic ray contamination after 
using the tagger described above. This is verified by in- 
specting distributions such as the timing of the COT drift 
chamber hits and the acolinearity of the lepton tracks. 
None of the candidate events are consistent with cosmic 
rays in any of the distributions. 

We calculate a cross section for each final state 
') using the formula 



N - N, 



bck 



a X JCeff' 

where N is the number of candidate events, N^ck is the 
expected number of background events, and a is the ac- 
ceptance. Assuming equal rates for the fJ,~^n~ and e+e^ 
processes, a combined cross section is found by forming 
a joint likelihood for the three final states, which is the 
product of the Poisson probabilities to observe N events 
in each final state. The joint likelihood is a function of 
the signal cross section, and all the systematic uncer- 
tainties are nuisance parameters. The likelihood is max- 
imized as a function of the signal cross section, and the 
uncertainty is obtained by multiplying this likelihood by 
a prior that is flat for positive cross sections and find- 
ing the shortest interval containing 68% of the integral. 
We treat the background and luminosity systematics as 
correlated and the acceptance systematics as uncorre- 
lated. The combined cross section for one lepton flavor 
is found to be a{pp Pjjp — > pi'^i^p) = 0.24^^ :[q pb 
for Mu > 40 GeV/c^ and \rii\ < 4, which IS m good 
agreement with the LPAIR prediction of 0.256 pb. 

Some of the kinematic properties of the candidate 
events are given in Table |T1 where p^{l) and p^{2) are 
the lepton transverse momenta, A(f)u is the difference in 



the azimuthal lepton angles (i.e. 180° minus A(j)u is the 
deviation from back-to-back in the transverse plane) and 
PT{£i) is the pT of the lepton pair. The resolution of the 
lepton transverse momenta is approximately 3.5(1.4)% 
for electrons (muons). All of the events have lepton pairs 
that are back- to-back in azimuth with low PTilt) values, 
which is expected for 77 t^t~ events. Figures 2(a) 



TABLE I; Properties of the eight exclusive dilepton events, 
in order of Mu- 



Final state 


Mu 




p^(2) 


180° - A<^« 






(GeV/c^) (GeV/c) 


(GeV/c) 


(degrees) 


(GeV/c 


e+e-(CC) 


40.9 


20.4 


20.1 


0.38 


0.26 


e+e-(CC) 


49.3 


24.5 


24.6 


0.37 


0.21 


e+e-(CP) 


50.4 


20.5 


20.2 


0.05 


0.31 


e+e-(CP) 


56.3 


24.8 


24.9 


0.48 


0.24 




58.6 


24.1 


24.4 


0.17 


0.32 




66.0 


31.8 


31.3 


0.75 


0.65 


e+e-(CP) 


67.1 


24.1 


24.0 


0.51 


0.24 


e+e-(CP) 


75.6 


34.1 


33.1 


0.23 


1.01 



and |2(b)] show the dilepton invariant mass and 180° mi- 
nus A(f>u distributions for the data together with the 
QED spectrum from lpair and the GEANT detector sim- 
ulation. A good agreement with the data is observed. 

No events pass our exclusive Z £^£^ selection crite- 
ria, therefore we place an upper limit on the cross section 
of exclusive Z production at the Tevatron. We sum the 
three final states to give ^ N = 0, N^ck — 0.66 ± 0.11, 
and ax //:e//xBR(£+£-) = 3.22 ± 0.38 pb^^. Here we 
have used BR(£+£~) = 3.37% as the branching fraction 
of the Z to decay to one lepton fiavor pair. We use a 
Bayesian limit technique to set an upper limit on the 
exclusive Z cross section of aexci{Z) < 0.96 pb at 95% 
confidence level. We also set an upper limit on the differ- 
ential cross section with respect to yz at yz—0 {^\y=o) 
using the theoretical prediction of the yz distribution fHI . 
We take 0.257 as the ratio of ^\y=o to crexci{Z) and find 

^|y=o < 0.25 pb at 95% confidence level. 

At hadron colliders the lepton kinematics in 77 
£'^£^ events determine the momenta of the for- 
ward (anti)protons through the relation ^(pi(2)) — 

75 Ei=i,2PTe^^"^'''' 0,11, where C(pi(2)) is the frac- 
tional momentum loss of the forward (backward) hadron. 
In principle this relation could be used to calibrate both 
the momentum scale and resolution of forward proton 
spectrometers. In our eight candidate events, only one - 
that with M^^ = 66.0 GeV/c^ - was from a period when 
the Roman pot spectrometer was operational and with 
^(p) in its acceptance; a track is observed, as expected 
for exclusive dilepton production. This is an encouraging 
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FIG. 2: (a) The dilepton invariant mass distribution, and (b) 
the distribution of 180° minus the difference in the azimuthal 
lepton angles for the data and the lpair prediction with the 
GEANT detector simulation, scaled to account for acceptance 
and luminosity. 



sign that if enough exclusive dilepton events are observed 
at the large hadron collider (LHC), they may be used to 
calibrate forward proton spectrometers 

In conclusion, we have observed exclusive production 
of high mass (Af« > 40 GeV/c^) e+e~ and ji" pairs 
and measured a cross section that agrees with QED ex- 
pectations. We observe no candidates for exclusive Z pro- 
duction and put an upper limit on the photoproduction 
of Z at a level k, 3, 000 times higher than SM predictions. 
It should be noted that an observation at the LHC, where 
the SM cross section is predicted to be 13 fb [5], is more 
promising. 
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